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1. Introduction 

Except for a few prominent counterexamples El Ej 5 the low energy physics 
of practically all chaotic quantum systems is governed by the universal spectral 
correlations of Wigner and Dyson's random matrix (RM) ensembles 0] . Yet in spite of 
its ubiquity, and notwithstanding a number of significant recent advances 00 El 
the correspondence above is not yet fully understood theoretically. Specifically, 
the 'non-perturbative' aspects of the problem — which manifest themselves, e.g., in 
the low energy profile of spectral correlations — are not under quantitative control. 
Some time ago, the introduction of a field theoretical approach, similar in spirit 
to the CT-models of disordered fermion systems, added a new perspective to the 
problem ^1 ^] . This so-called 'ballistic cx-model' describes chaotic systems in terms 
of a field theory in classical phase space. Remarkably, it provides a faithful description 
of RM spectral correlations already on the most elementary mean field level where 
fluctuations inhomogeneous in phase space are neglected; 'all' that remains to prove 
universality is to show that these inhomogeneities indeed become inessential in the 
long time limit — an expectation backed up by the long time ergodicity of chaotic 
systems. 

Unfortunately, however, this latter task soon proved to be excruciatingly difficult. 
In this paper we shall concentrate on the perhaps most serious of these problems, 
the seeming incapability of the new approach to correctly describe even the lowest 
order quantum interference corrections ('weak localization corrections' in the jargon 
of mesoscopic physics) to physical observables: In a semiclassical manner of speaking, 
'quantum interference' is a process wherein two initially identical — modulo the 
notorious uncertainty introduced by the non-vanishing of Planck's constant - 
Feynman trajectories split up and later recombine to an overall phase coherent 
structure (see figure ^) . This mechanism is at the root of practically all quantum 
phenomena distinguishing disordered or chaotic quantum systems from their classical 
limits. It is closely tied to the notion of the Ehrenfest time — the time it takes 
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for the separation of two trajectories to grow from Planck scales to macroscopic 
scales. Irritatingly, however, the field theory formalism appeared to be incapable of 
describing the initial ft-uncertainty triggering these phenomena. Deferring for a more 
detailed discussion to section [3] below, let us try to outline the essence of the problem: 
loosely speaking, the field degrees of freedom of the ballistic er-model describe the joint 
propagation of retarded and advanced Feynman amplitudes along classical trajectories 
in phase space. Previous works effectively did not allow for deviations between the two 
amplitudes. At this level of approximation, the retarded and the advanced reference 
path are strictly identical and the fi-quantum uncertainty essential to initiate the 
formation of quantum interference corrections is absent. Equally important, points 
in phase space belonging to different classical trajectories remain uncorrelated. This 
implies that the theory will not be able to describe the relaxation into the uniform 
mean field configuration (i.e. will not be able to predict RMT behaviour.) 

A phenomenological solution to this problem was proposed by Aleiner and 
Larkin (141 1151 116) . Building on the insight gained in previous work, they added 
a diffusive contribution (formally, a second order elliptic operator) to the action 
of the model. Multiplied by a coupling constant of 0(h), this term introduced 
a sufficient amount of 'fuzziness' to the problem to initiate quantum interference 
processes. Although the extra contribution to the action could not be derived from 
first principles, AL argued that it ought to be present on physical grounds (viz. to 
mimic the diffractive aspects of the propagation of quantum states.) 

It is the purpose of this paper to demonstrate that, in fact, no diffraction terms are 
needed to describe quantum interference within the framework of the ballistic <7-model. 
Our analysis will hinge on an aspect of the theory that has been noticed before ^2] 
yet did not receive sufficient attention: the cr-model is not a local field theory in 
phase space; by construction, and in accord with the principles of the uncertainty 
relation its maximal resolution is limited to Planck cells of extension ~ tV, where / 
denotes the number of degrees of freedom. We will show that this non-locality suffices 
to describe quantum interference in far-reaching analogy with recent semiclassical 
approaches [SJ m to the problem. t In recent work [T^J, similar ideas have been 
applied to compute (in a non field theoretical setting) weak localization corrections 
of a quantum map (viz. the standard map or kicked rotor). 

Specifically, we will consider the spectral two-point correlation function R2 (li) at 
energies u> larger than the single particle spacing A. We will show that the expansion 
of i?2 in the small parameter = (ttlo/ 'A) -1 agrees with the prediction of RMT. 
(In a manner that largely parallels our present analysis, the same result has recently 
been obtained by periodic orbit theory 22-) The extensibility of the analysis to the 
perturbatively inaccessible regime s < 1 remains an open issue. 

The rest of the paper is organized as follows: To facilitate the comparison with 
the field theoretical formalism, we begin by reviewing some of the recent developments 
in the semiclassical approach to quantum chaos (section^. In section[3]we turn to the 
field theoretical approach and apply it to the perturbative expansion of the two-point 
correlation function. We conclude in section 

f It is due to mention, though, that our analysis, too, necessitates the ad— hoc addition of an extra 
contribution to the action of the native model. Yet, in a sense to be qualified below, this term serves 
purely regulatory purposes. Coupled to the theory at a strength parametrically weaker than that of 
the AL term, it does not affect the dynamics at times t < tj> 
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2. Semiclassical Background 

We are interested in the behaviour of globally hyperbolic (chaotic) quantum systems 
at time scales t larger than the ergodic time £ C rgt yet smaller than the Heisenberg 
time <h = 2nh/A. (The first condition implies that non-universal aspects of the 
classical dynamics are inessential, the second that concepts of perturbation theory (in 
the parameter r = t/tn) are applicable.) 

To describe correlations in the spectrum of the system, we consider the two-point 
correlation function 

R 2 {uj)=A 2 (p{E + uj/2)p(E-uj/2)) E -l (2.1) 

and its Fourier transform 

K(t) = ^y^e-^^M, (2.2) 

the spectral form factor. Here, p(E) is the energy dependent density of states (DoS) 
and (. . .) e denotes averaging over a sufficiently large portion of the spectrum centered 
around some reference energy Eq. 

In semiclassics, the spectral form factor is expressed as 

77' 

where ^ s a double sum over periodic orbits 7 and 7', S 7 the classical action of 

the orbit, T 7 its revolution time, and Ay its classical stability amplitude. 

Before turning to a more detailed discussion, let us briefly summarize the main 
results recently obtained for the semiclassical form factor: For times r < 1, K sc can 
be expanded in a series in r. As shown by Berry |19|. the dominant contribution to 
this expansion Ksc = 2r, is provided by pairs of identical (7 = 7') or mutually time 
reversed (7 = T7') paths. (Throughout we focus on the case of time reversal and spin 
rotation invariant systems — orthogonal symmetry.) 




Figure 1. Cartoon of a pair of topologically distinct paths, (7,7') contributing 
to the first quantum correction to the spectral form factor. Notice that 7 and 
7' differ in exactly one intersection point (crossing vs. avoided crossing). Inset: 
blow-up of the intersection region. 

f Formally, t crg is defined as the inverse of the first non-vanishing Pcrron-Frobcnius eigenmode. 
In fact, all our results can be generalized to general mixing rather than just uniformly hyperbolic 
systems. The point is that mixing implies ergodicity and non-integrability, and hence any mixing 
system will appear to have a constant global Lyapunov exponent when evaluated on time scales 

t ^> terg ■ 
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All corrections to the leading contribution hinge on the mechanism of 

quantum interference alluded to in the introduction. E.g., the sub-dominant 

(2) 

contribution, JQc , to the form factor is provided by pairs (7,7') that are nearly 
identical except for one 'encounter region'rj In this region one of the path self- 
intersects while its partner just so avoids the intersection (cf. figure^). (Alternatively, 
one may think of two trajectories that start out nearly identical, then split up and later 
recombine to form an interfering Feynman amplitude pair.) The two paths are, thus, 
topologically distinct yet may carry almost identical classical action ^1] . Specifically, 
Sieber and Richter [5] have shown that for sufficiently shallow self intersections 
(crossing angle in configuration space of 0(H)) the action difference \S 7 — 5 7 < < h. 
For these angles, the duration of the encounter process is of the order of the Ehrenfest 
time ts = A -1 ln(c 2 /?i), where A is the phase space average of the dominant Lyapunov 
exponent of the system and c a classical reference scale (see below) whose detailed value 
is of secondary importance. This identifies ts as the minimal time required to form 
quantum corrections to the form factor (as well as to other physical observables |14)V 
Throughout we shall assume t crg < £e < t < tji, where the condition t erg < 
is imposed to guarantee that for time scales t > t^, the system indeed behaves 
universally. (For t CIg > the time window £e < t < t elg is characterized by the 
prevalence of correlations that are non-universal yet quantum mechanical in nature.) 



Figure 2. Cartoon of three classes of orbit pairs that contribute to the expansion 
of the form factor at order r 3 . (The triple-encounter region shown in the two 
figures on the right is the analog of the Hikami hexagon familiar from the impurity 
diagram approach to disordered systems.) The existence of the middle pair docs 
not rely on time reversal invariance. 

Summation over all Sieber-Richter pairs leads to the universal result K sc ~ 
KsP + Ksc = 2r — 2t 2 , which is consistent with the short time expansion of the 
random matrix form factor 

K rm (t) °-=- 1 2r-rln(l + 2r). (2.3) 

At higher orders in the r-expansion, orbit pairs of more complex topology enter 
the stage. (For some families of pairs contributing to the next leading correction, 
see figure |2J) The summation over all these pairs — feasible under the 
presumed condition t CIg < r — obtains an infinite r-series which equals the series 
expansion of the RMT result p.HJl .S It is also noteworthy that both the topology 
of the contributing orbit pairs and the combinatorial aspects of the summation are 
in one-to-one correspondence to the impurity-diagram expansion |21j of the spectral 
correlation function of disordered quantum systems. 

I Notice that a path of duration t 3> terg generally contains many self intersections in configuration 
space . 

§ However, as is indicated by the notorious non— analyticity of Krm(t) at r = 1 1201 . the form factor 
at times r > 1 appears to be beyond the reach of semiclassical summation schemes. 
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Figure 3. The structure of the encounter region. The picture on the right 
shows how the parallelogram spanned by the four points evolves in time t\ , while 
its symplectic area us is conserved. 



Central to our comparison of semiclassics and field theory below will be the 
understanding of the encounter regions where formerly pairwise aligned orbit stretches 
reorganize. The analysis of these objects is greatly facilitated by switching from the 
configuration space representation originally used in to one in phase space 00 Ell- 
in the following we briefly discuss the phase space structure of the regions where 
periodic orbits meet. In section 13.41 we will compare these structures to the - 
somewhat different — field theoretical variant of encounter processes. 

(2) 

Considering the correction iv sc as an example, we note that the encounter region 
contains four orbit stretches in close proximity to each other (cf. figures □ EJ: two 
segments x(£i) and x'(ii) of the orbits 7 and 7' traversing the encounter region and 
the time reversed\\ x(ii + £2) and x'(ii + £2) of the trajectories reentering after one 
of the loops adjacent to the encounter region has been traversed (£2 is the duration 
of the loop traversal and t\ parameterizes the time during which the encounter region 
is passed). To describe the dynamics of these trajectory segments, it is convenient 
to introduce a Poincare surface of section S transverse to the trajectory x(ii). For 
a system with two degrees of freedom (a billiard, say), S is a two-dimensional plane 
slicing through the three-dimensional subspace of constant energy in phase space. We 
chose the origin of S such that it coincides with x(ii ) . Introducing coordinate vectors 
e u and e s along the stable and unstable direction in S, the three points x(ti + t%), 
x'(ti) and x'(ti +ta) are then represented by the coordinate pairs (u, s),(u, 0) and 
(0, s), respectively. (Notice that the trajectory x'/x' traverses the encounter region 
on the unstable (s = 0)/stable [u = 0) manifold thus deviating from/ approaching the 
reference orbit x.) 

The above coordinate system is optimally adjusted to a description of the two 
main characteristics of the encounter region: its duration i enc and the action difference 
Indeed, it is straightforward to show that the total action difference is simply 
given by the area of the parallelogram spanned by the four reference points in phase 
space, S y — Sy = us 0. As for the encounter duration, let us assume that the 
distance between the orbit points may grow up to a value c before they leave what 
we call the 'encounter region'. (It is natural to identify c with the typical phase space 
scale up to which the dynamics can be linearized around x(fi), however, any other 
classical scale will be just as good.) After the trajectory x has entered the encounter 
region, it takes a time t ln ~ A _1 ln(c/s) to reach the surface of section and then a 
time i out ~ A ln(c/tt) to continue to the end of the encounter region. (Here, A is 

|| In a standard position-momentum representation x = (q, p), time reversal is defined as x = 

(q, -p)- 
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the Lyapunov exponent of the system. Thanks to the assumption i erg <C t-g, A may 
be assumed to be a 'self averaging quantity', constant in phase space.) The total 
duration of the passage is thus given by t enc (u, s) = t out + t m ~ A -1 ln(c 2 /(us)). The 
action difference of orbit pairs contributing significantly to the double sum must be 
small, |5 7 — S 7 '| = us < h. Consequently, t cnc > t~E = \~ l \n.{c 2 /h), where t~E is 
the Ehrenfest time introduced above. (Notice that both S* 7 — Sy and t enc depend 
only on the product us. While the individual coordinates u and s depend on the 
positioning of the surface of section, their product us is a canonical invariant and, 
therefore, independent of the choice of S.) 

Having discussed the microscopic structure of the encounter region, we next need 
to ask a question of statistical nature: given a long periodic orbit 7 of total time t, 
what is the number N(u, s,t)duds of encounter regions with Poincare parameters in 
[u,u + du] x [s,s + ds}1 (To each of these encounter regions there will be exactly 
one topologically distinct partner orbit 7' that is identical to 7 in all other (N — 1) 
encounters. Thus, N(u, s,t)duds is the number of Sieber-Richter pairs for a given 
parameter configuration and J dudsN(u, s, t) is the total number of Sieber-Richter 
pairs.) Since the times t\ and t 2 defining the two traversals of the encounter region 
are arbitrary (except for the obvious condition \t x — t 2 \ > t cnc ), N is proportional 
to the double integral N(u, s,t)duds oc 5 J" u 2 _ tl i >t dt\dti P rct (u, s,ti)duds. The 
integrand, P rct is the probability to propagate from the point (0, 0) in the Poincare 
section to the time-reverse of (u,s) in time £2- Since ti > £e > ^crg, this probability 
is constant and equals the inverse of the volume £1 = 2nhtii of the energy shell, 
P ret (u, s,t2) = fl^ 1 . Thanks to the constancy of P re t, the temporal integrals can 
be performed and we obtain N(u,s,t) oc t(t — 2i cnc ) /2fi. The normalization of N 
is fixed by noting that the temporal double integral weighs each encounter event 
with a factor t cnc . The appropriately normalized number of encounters thus reads 
N(u, s) = *^2t 2t °n ^ • Substitution of iV(it, s, t) into the Gutzwiller sum obtains 

K(2) ( T ) =Y,\ A i\ 2s ( t - J 1 ) f dudsN{u,s,t)2cos(us/H) 

= — — / duds I ; r — 2 ) costus/h) l = — 2t 2 . 

2nhJ_ c \t cnc (u,s) J v 1 1 

where we used the sum rule S 7 l^7l 2 ^( T — f~) = t of Hannay and Ozorio dc 
Almeida [22] and noted that in the semiclassical limit the first term in the integrand 
does not contribute (due to the singular dependence of t cnc on H.) 

Before leaving this section, let us discuss one last point related to the semiclassical 
approach: the analysis above hinges on the ansatz made for the classical transition 
probability Pt(x,x') between different points in phase space. Specifically, a naive 
interpretation of ergodicity — P t (x, x') = SI -1 = const, for times t > t eTg — is too 
crude to obtain a physically meaningful picture of weak localization. One rather has 
to take into account that the unstable coordinate, u(t), separating two initially close 
(w(0) <C c) points x and x' grows as u(t) ~ u(0) exp(At). For sufficiently small initial 
separation, the time it takes before the region of local linearizability is left, 

|t E (x,xO = \ln^ y (2.4) 

may well be larger than i erg . This is important because during the process of 
exponential divergence, the probability to propagate from x to the time-reverse x' 
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is identically zero. (Simply because the proximity of x and x' implies that x and x' 
are far away from each other.) Only after the domain of linearizable dynamics has 
been left, this quantity becomes finite and, in fact, constant: 

P t (x,x') = ^G(t-t E (x,x')). (2.5) 

This concludes our brief survey of the semiclassical approach to quantum coherence. 
We next turn to the discussion of the field theoretical formulation and its structural 
parallels to the formalism above. 

3. Field Theoretical Formulation 

3.1. Definition of the Model 

The ballistic cr-model is defined by a functional integral Z{lo) = J VTe~ s ^ extending 
over field configurations T(x) in classical phase space. Its action is given by S[T] = 
S [T] + S Ieg [T], where 

So[T] = ihj (dx) t^^-'AT + T-^T]) (3.!) 

is the action of the 'native' model |13| and S reg a regulatory contribution to be 
discussed momentarily. Here, f(dx) = (2nK)~f +1 J dx5(Eo — H{x)) is the integral 
over the (2/ — 1) -dimensional shell of constant energy Eq\, T = T(x), cf r and A 
are matrices whose internal structure will be discussed momentarily, H the classical 
Hamilton function of the system, and lo the scale at which we are probing the spectrum. 
(Within the field theoretical approach it is preferable to work in energy rather than 
in time space.) Importantly, all products appearing in the action l|3.1fl have to be 
understood as Moyal products, 

(/ 5 )(x)^e^ 9 - ra -/( Xl ) 3 (x 2 )| xi=x2=x , 

where the matrix / is defined through x T ix' = q ■ p' — p ■ q'. For later reference we 
note that the Moyal product affords the alternative representation 

(/5)(x) -/ ww efxr/X2/(x+xi) - 9(x+X2) - (3 - 2) 

Equation (|3.2|) makes the 'non-locality' inherent to the action of the ballistic cr-model 
manifest: all products involve a coordinate averaging over Planck cells of volume 
~ fv . As we shall see below, this non-locality encapsulates essential aspects of the 
semiclassical dynamics discussed in the previous section. 
The second contribution to the action 

SregPI = y (dx) tr{d l T(x)d l T- 1 (x)) (3.3) 

serves to damp out singular field configurations (Unlike with most other field theories, 
the action of the unregularized model, governed by the generator [H, ] of unitary 
quantum dynamics, does not have the capacity to self-regularize.) In | Appendix A| we 
will argue that a coupling constant g TeE ~ f? suffices to stabilize the theory. Coupled 
to the theory at this strength, the action S lcg does not yet influence the dynamics on 
the physically relevant times £e- This stands in contrast to the theory of AL where a 

f See |Appcndix A] for details on the definition of this integral. 
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second order derivative term (similar in structure to (|3.3|) but with coupling constant 
<?reg ~ K) actively governed the dynamics at times t ~ t-g. 

In the original references, the ballistic cr-model was introduced as a 
supersymmetric field theory. However, for the purposes of our present analysis it will 
be more convenient to employ the simpler formalism of the replica trick. Within this 
approach, the matrices T(x) S Sp(4R)/Sp(2i?) x Sp(2i?) act in a tensor product of R- 
dimensional replica space, a two-dimensional space distinguishing between advanced 
and retarded propagators (ar-space) and a two-dimensional space (tr-space) whose 
presence is required to account for the time reversal invariance of the system Here, 
Sp(2i?) is the 2i?-dimensional symplectic group and A = erf® l tr (g>lfl, where 1_r is the 
i?-dimensional unit matrix. While the use of replicas bars us from performing non- 
perturbative calculations, it significantly facilitates the comparison to the semiclassical 
analysis above. 



3. 2. Two-Level Correlation function 

Our goal is to compute the two-level correlation function i?2(w). Expressed in terms 
of single particle Green functions G^ = (E ± iO — -ff) -1 , 

R 2 {u) = |^Re <tr(G + (£ + w/2)) tr(<T(£ - u/2))) E c , 

where (. . .}_e. c denotes the (connected) average over an energy interval much larger 
than the inverse of the smallest time scales in the problem (the inverse of the Lyapunov 
time, say.) Within the replica formalism, R2 is obtained by a two-fold differentiation 
of the partition function w.r.t. the energy parameter:^ 

M*) = -j £m ^Red 2 s Z(s), 

where the dimensionless variable s = ttuj/A. As long as we restrict ourselves to 
perturbative operations, i.e. an expansion of the two-level correlation function in a 
series 

00 

i? 2 ( S )^ 1 Re^c n (i S +)-", (3.4) 

71=2 

the replica limit R — ► is well-defined. A straightforward Fourier transformation, 
K(t) = 7T _1 J ds e _2lST i?2(s)) shows that the coefficients c„ are related to the 
coefficients d n of the spectral form factor K (r) = dnT n through 

(-2)™ 

d n = j C n +1- (3.5) 

IV. 

In fact, however, there are much more far-reaching analogies between the temporal 
and the frequency representation of spectral correlations: at every given order 
n various topologically distinct families of orbit/partner orbit pairs ('diagrams') 
contribute to the coefficient d n . Likewise, the expansion coefficients c n obtain as 

I This follows from the fact that (by construction) 

Z{ui -0) 3 ) = <dct[iG+(£; + wi)] B dct[iG-(£; + a)2)] fl >i5. 
Using that In 2 = lirnR_,o( zi? — 1)/-R, it is then straightforward to verify that 

Urn •^■Re92 1 _ Wa £(wi -0)2) = -Re(tr G+(£ + wi)tr G~ (E + oj 2 )) e ,c = -2(tt / A) 2 fl 3 (wi -wa). 
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sums of Wick contractions of the generating functional Z(ui). We shall see that 
there is an exact correspondence between field theoretical and semiclassical diagrams 
(both in topological structure and numerical value) which simply means that the 
two approaches describe spectral correlations in terms of the same semiclassical 
interference processes. 

3.3. Quadratic Action 

We next turn to the actual expansion of the field integral. For this purpose, we shall 
employ the so-called 'rational parameterization' of the coset-valued field T. This 
parameterization is defined by T = 1 + W, where 

w ={b M 

\ / ar 

is a matrix that anti-commutes with the matrix A introduced above. Its off-diagonal 
blocks take values in the Lie algebra of Sp(2i?), i.e. they satisfy the constraint 

= B T = (icrf ® 1r)B t (ia^ ® 1r) _1 - The principal advantage of the rational 
parameterization is that the Jacobian associated to the transformation from the T- 
matrices to the linear space of i3-matrices is unity: J T>T = J T>B. 

Substituting this representation into the action, we obtain a series S[B] = 
S^ 2n ^[B], where S^ 2 ™) is of 2n-th order in B. Let us begin by discussing the 
unregularized quadratic action 

S » )[B] = ~Th\ {dx) tr [jBt (w+ ~ [H > ]) s] ' 

where [H, ] is the generator of quantum time evolution. Very little can be said about 

(2) 

this generator in concrete terms which means that the action Sq does not qualify 
as a 'reference point' of a perturbative expansion scheme. (Indeed, notice that the 
projection |a)(a| onto an eigenstate \a) of the Hamilton operator H is annihilated by 
[H, ]. This means that the quantum evolution operator possesses a large number 
of (nearly) unstable 'zero modes' whose action is damped only by the frequency 
parameter to.) 

Much more is known about the generator {H, } of classical dynamics, where 
{/, g}(x) = (?x 1 -fd X2 /(xi)<7(x2)| _ x _ x is the Poisson bracket. Expanding the Moyal 
commutator, 

[H, B](x) - ih{H, 5}(x) + 0((hd x ) 3 B( x )), 

we notice that the quantum generator [H, ] differs from its classical counterpart {H, } 
by the presence of higher order derivative terms. In |Appendix A| it is shown that the 
quadratic regulatory action 

S%l[B}=g Ieg [ (dx) tr(9^t(x)a 4 B(x)) 

suffices to damp out higher derivatives and hence effectively projects the quadratic 
action onto its classical limit. Assuming regularization in the above sense, our further 
discussion will be built on the action 

Si? ^ = \j ( dx ) tr [S t (x)(^S)(x)] , (3.7) 

where £ w = —iu/h — {H, }. Throughout, the operator = £~r eg will play 
an important role. Here, the subscript 'reg' indicates that acts in a space of 
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functions coarse grained over cells in phase space of 'volume' (h 2 /a)f, where a is 
some classical reference scale of dimensionality 'action' whose detailed value will not 
be of much concern. Importantly, P w is not strictly inverse to the bare Liouville 
operator (i.e. the Liouville operator acting in the space of unregularized functions.), 
x') ^ <5(x — x'). Rather, the time Fourier transform, P t (x, x') = <5(x — x(f)) 
can resolve the definite dynamical evolution generated by the Liouville operator only 
up to time scales 

2 

t E = A^ln^ ~ 2tE. 

AT 

Thereafter, the uncertainty in the resolution of the boundary conditions (the effect of 
smoothening) renders the dynamics unpredictable, i.e. 




the crossover between the two regimes takes place over time scales ~ max{A?E, ^crg}, 
where A^e <C £e is the uncertainty in £e caused by an eventual non-uniformity of the 
Lyapunov expansion. § (Notice that in previous discussions of the ballistic tr— model, 
the propagator P was mostly identified with the Perron-Frobenius operator, i.e. an 

object that describes relaxation into a uniform configuration, -Pt(x, x') = rg const, 
over classically short times. However, it is not at all clear how this behaviour 

may be reconciled with the indispensable condition that P t (x, x') t<tE i*' x ** Q: for 
|x — x'| = 0(h), the propagator must be able to resolve fine structures in phase space 
over times ~ tg parametrically larger than the relaxation time of the Perron-Frobenius 
operator. In contrast, equation ljii.81) is motivated by the structure of the action, and 
does resolve the classical phase space dynamics up to the Ehrenfest time. In fact, we 
will see that the scale ?e > £e does not explicitly enter the results of the theory. The 
reason is that, due to a conspiracy of quantum non-locality and chaotic instability, 
the dynamics becomes effectively irreversible instantly after the Ehrenfest time (on 
time scales of the order of the inverse Lyapunov exponent). Thus, at a time where the 
artificially introduced smearing would become virulent, the theory has long become 
quantum-unpredictable.) 

S-4- Perturbative Expansion 

We now have everything in store to proceed to the perturbative expansion of 
the functional integral. The dominant contribution to the series (|3.4f) obtains by 

§ The results above apply to uniformly hyperbolic systems. In the case of non-uniform hyperbolic 
systems, local fluctuations in the Lyapunov expansion rate A(x) need to be taken into account. The 
logarithmic mismatch j/(x, i) = \a(u(t) / u{0)) between two trajectories starting at x and x + «(0)e u , 
respectively, grows as y = A(x(t)). (e u is the locally most unstable direction in phase space.) Due to 
inhomogeneities in the expansion rate, y(x, t) is a fluctuating quantity with mean y(t) and a certain 
width Ay(t). Importantly, an upper bound on fluctuations in y is imposed by Oseledec's theorem 1251 
which states that the phase space average A of the Lyapunov expansion rate equals the long— time 
expansion rate of individual trajectories almost everywhere: y(x,t)/t — > A for t — > oo for almost all 
x. Consequently, Aj/(i) ~ t a grows at a rate a < 1. (E.g., the model of statistically independent 
Gaussian fluctuations of the local expansion rate employed by AL 1141 leads to a = i.) By definition 
of £e> an phase space distribution of initial extension ~ h 2 f has expanded to classical dimensions 
when y(t) = Afg. Defining Atg as the time uncertainty in tj; (due to fluctuations in the local 
expansion rate), we obtain the estimate A£e ~ Aj/(tE)/A ~ tS. This means that Ats/iE ~ t a ~ 1 
vanishes in the semiclassical limit. For finite h, the effective relaxation rate of the system is set by 
max{At E , tag}- 
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integration over the quadratic action: 



i Re Iim ±-d 2 s {detP u ) 2R2 = Red 2 s lndet^J 1 ) Re(i S +)- 2 . 



(3.9) 



This result implies (cf. equations (|3.4I I3.5( l) di — 2 in accord with the semiclassical 
analysis. || 

To compute higher order terms in the expansion we need to consider the non- 
linear contributions S^ 2n>2 ^ to the action. Substituting the representation (|3.6|l into 
the action l|3.1(l we obtain 

S {2n) [B] = \ J (dx)tr [(-B t B)"- 1 S t £ tJ B] . (3.10) 

By elementary power counting, each matrix B scales as (symbolic notation) ~ 
(£ w )~ 5 uj-2 ^ s~2. This implies that each vertex S^ 2 ™) contributes a factor 
~ (B^ B)^ 1 B^ C U B ~ s~ n+1 to the functional integral. Specifically, the dominant 
correction (~ s -3 ) to the leading contribution l|3.9|) obtains by first order expansion 
in the vertex 



R ( i\s) = - Re (dx)<tr(Bt B Bt jCwjB )) B . 



(3.11) 



We emphasize again that all products of B-matrices have to be understood as Moyal 
products. To obtain a convenient representation of the product of more than two 
of these matrices, we iteratively apply the prototype formula equation 13.211 . A 
straightforward calculation then yields the general product formula 

(A,... A 2n )(x) = /I] T^yj ^^—^(x + Xl ) . . . A 2 „(x + x 2n ), 

i— 1 ^ 1 

where the bilinear form S(xi, ■ ■ -x 2 „) = 2 Yli<j '( — ) l+:,x f ^ x j- Below, we will apply 
this formula to the fields B of the theory. In |Appendix A| we show that in this 
case, all energy coordinates Ei get locked. Here, we assume a coordinate choice 
x = (£?, t, y) where E is an energy variable, r its canonically conjugate time-like 
variable (a coordinate parameterizing the Hamiltonian flow through x) and y a 
(2/ — 2) -dimensional vector of coordinates transverse to the flow. Further, fluctuations 
in the time-like variables Tj are negligible. Introducing the shorthand notation 
/ = (nh)~f +1 f dxi 5(Ei — E )S(Ti), we thus have 

(B^B) 2n ( X )= f e^ Xl <- X2 ") J Bt(x + x 1 )...B(x + x 2 „). ( 3 . 12 ) 

J xi,...,x 2 „ 

|| It is worthwhile to notice that the agreement between semiclassics and field theory does not pertain 
to times t < t CTg : for these times, short periodic orbits traversed more than once influence the 
behaviour of the form factor. For reasons that only partly understood, the cr-model fails to correctly 
count the integer statistical weight associated to the repetitive traversal of periodic orbits. The 
essence of the problem 1201 is that the degrees of freedom of the a— model (the B's) describe the 
joint propagation of amplitudes locally paired in phase space. However, an n-fold repetitive process 
is governed by the local correlation of 2n Feynman amplitudes. Perturbative approaches to the 
problem fail to correctly describe these correlations. Interestingly, a non-perturbative evaluation of 
the functional integral — feasible in the artificial case of periodic orbits with unit monodromy matrix 
— leads to the correct result (M. R. Zirnbauer, unpublished). 
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Using this representation in l|3.11fl . applying the contraction rules (|lj.l|) discussed 
in |Appendix B| and taking the replica limit we obtain 



(s) = Re^Y d 2 s J (dx) J e^ Xl — x ^P w (x + x l5 x + x 3 )£ w , X4 P,(x + x 4 , x + x 2 ), 



4 3) 

/Xl,...,X4 

where the coordinate subscript in £ WiX indicates the argument on which the Liouvillian 
acts. The physical meaning of this expression is best revealed by switching to the 
Fourier conjugate picture. Inserting the definition l|2.2(l of the form factor, we obtain 

12 



K^(t) = -2t 2 — /(dx) / e i s ^<-<^ 
ta J ix lr .,x 4 



/ dt'P t - t > (x + xi, x + x 3 )£ ( / iX4 P t / (x + x 4 , x + x 2 ), 
Jo 



(3.13) 



where Ct = dt — {H, }. equation H3.13|l makes the analogies (as well as a number 
of differences) between the semiclassical and the field theoretical description of 
quantum corrections explicit: central to both approaches are two semi-loops shown 
schematically in figure ^ I n either case, the proximity of these loops is controlled 
by phase factors which contain the coordinates of the end points (in a canonically 
invariant manner) as their arguments. However, unlike with semiclassics, equation 
<|3.13[1 does not relate the unification of the two semi-loops to specific periodic orbits. 
Rather, the two halves are treated as independent entities, each described in terms of 
its own probability factor P. Relatedly, the phase factor controlling the proximity of 
the terminal points does not correspond to the action difference between two orbits. 

The result obtained for K^(t) in equation (|3.13(1 critically depends on the 
behaviour of the propagator Pt at times t ~ t E (x, x') ~ t E , cf. equations (|2.41 12. 5| . 
Specifically, we shall use that d(Pt(x, x') = f2 _1 (5(i — t E (x,x.')), where S(t) is some 
smeared (5-function whose detailed functional structure is not of much importance. 
(All we shall rely upon is J dt5(t)f(t) ~ f(t) for functions that vary slowly on the 
scales where 5(t) varies.) We also note that the Poisson bracket {H, /}(x) ~ <9 t /(x) 
effectively differentiates along the trajectory through x. However, the time t E (yi, x ') = 
t E (y,y') defined in equation l|2.4|l depends only on the coordinates transverse to the 
trajectory. This implies {H, Pt(x, x')} = {H, P t (t E (x, x'))} — 0. We thus conclude 
that the action of Ct on the function P is given by £ Xj iPf(x, x') ~ 1 <5(t — t E (x, x')) ~ 
Q,~ 1 S(t — tE)- To understand the meaning of the second approximation, notice that it 
takes a time tE before the bulk of the Planck cell to which the points x and x' belong 
has grown to classical scales c. For times t > tE, the fraction of the Planck cell which 
has not yet acquired macroscopic dimensions shrinks exponentially on the classical 
Lyapunov scale A -1 ^ tE- This means that ts(x, x') ~ up to an insignificant 
uncertainty of 0(X^ 1 ). Using these results, as well as the normalization relations 
/ (dx) = % and f xa e iS(xi,..., X2 „) _ ^ we bt am 



dt'Q{t - t' - t E )8(t' - t E ) 



K&(t) ~ -2r 2 — f (dx) f e iS(*i,...,**) I 

tft J ix I ,...,X 4 JO 

= -2T 2 e(^-2^E) 

in agreement with the result of the semiclassical analysis. 
3. 5. Higher Orders of Perturbation Theory 

What happens at higher orders in perturbation theory in the parameter s -1 ? Before 
turning to the problem in full, it is instructive to have a look at the zero mode 
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approximation to the model. The action of the zero mode configuration — formally 
obtained by setting T(x) = T = const. — is given by S [Q] = is + tr(af : Q)/4:, where 
we have used the standard |27| notation Q = T -1 AT. Parameterizing the matrix 
T = 1 + W as in 1|3.6|) , an expansion in the generators B obtains the expression 

oo 

S [B] = S 2n) i B l S^ n) [B] = -is+tr(-B t B)™. (3.14) 

n=l 

It is known |27) that, term by term in an expansion in a , the zero mode functional 
reproduces the RMT approximation to the correlation function i?2(s). Second, there 
exists a far-reaching structural connection between the perturbative expansion of the 
zero mode theory on the one hand and the Gutzwiller double sum on the other. (In 
fact, the correspondence Gutzwiller sum <-> zero dimensional c-model «-» RMT played 
a pivotal role in the proof that the semiclassical expansion coincides with the RMT 
result. [Til) 

More specifically, to each term contributing to the Wick contraction of 

((sW[B]r*(S^[B}r°...) (3.15) 

there corresponds precisely one semiclassical orbit /partner orbit pair (or 'diagram'). 
By power counting, this diagram contributes to the correlation function at order 
s~ 2 ~^" m "( n-1 ' ) . For every value of n = 2,3,..., it contains m n encounter regions 
where n orbit segments meet and ^ n nm n inter-encounter orbit stretches. The 
topology of the diagram is fixed by the way in which the B matrices are contracted. 
(For example, the first of the diagrams shown in figure[2]corresponds to the contraction 
(1-3,2-6,4-8,5-7) of tr(B^BB^B)ti(B^BB^B), the second diagram to the 
contraction (1 - 4, 2 - 5, 3 - 6) of tr(B^ BB^ BB^ B), etc.) Importantly, the minimum 
time required for the buildup of a diagram (i.e., the time required to traverse the 
J2 n TO ™ encounter regions) is given by £e x J2 n nm n- 

Turning back to the full problem, let us consider the analog of the zero dimensional 
expression <|3.15[) . 

((SW[B]) m * {S (6) [B}) m3 ...), (3.16) 

where S^ 2n ^ is given by (f 3 . 1 Cj|) and the average is over the full quadratic action (j3.7|l . 
It is natural to expect that the unique correspondence between Wick contractions and 
semiclassical diagrams carries over to the full model. If so, individual contractions 
should vanish/reduce to the universal RMT result for times shorter/much larger than 
t < £e x Yin nm n- In section l3~4l this correspondence was exemplified for the simplest 
non-trivial example, the Sieber-Richter diagram (S^[B}). 



Figure 4. Two representatives of the 'clover leaf diagram class contributing to 
the form factor at 0(t 3 ). Discussion, sec text. 

Perhaps unexpectedly, the straightforward one-to-one correspondence outlined 
above does not pertain to higher orders in perturbation theory. To anticipate our main 
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findings, it turns out that at order (s~ 4 <-> t 3 ) in the series expansion, propagators 
of short duration P«t E — absent in the (t 2 <-> s~ 3 ) term considered above — begin 
to play a role. This implies that individual contractions may relate to more than one 
semiclassical diagram class. Nonetheless, integration over all time parameters obtains 
a universal result. 

By way of example, let us consider the (1 — 3, 2 — 6, 4 — 8, 5 — 7) contraction of 
(tr(B t BB^ B)tr(B^ BB^ B)) . For generic values (U ~ t H > *e) of the time arguments 
carried by the four resulting propagators the contraction corresponds to the orbit 
pair shown in figure [21 left. However, the integration over times U also extends over 
exceptional values where one of the two propagators connecting the two encounter 
regions ((2 — 6) or (4 — 8)) is of short duration < t^- Such a short time propagator 
connects two distinct vertices.^ This results in a structure as shown in figure 01 
right, where the two clusters of dots indicate the eight phase space arguments of 
the -B-fields, the straight line-pair represents the short propagator, and the box 
indicates that all phase space points lie in a single encounter region. Evidently, this 
structure corresponds to a pair of orbits visiting a single encounter region twice. 
Diagrams of this structure are canonically obtained by contraction of a 'Hikami 
hexagon' tr(B* BB^ BB^ B), as indicated in figure^left. Fortunately, the absence of a 
unique assignment to semiclassical orbit families, does not significantly complicate the 
actual computation of the diagrams: closer inspection shows that taking the Liouvillc 
operators involved in the definition of the Hikami boxes into account and integrating 
by parts, we again obtain the universal zero-mode result. 

Summarizing, we have seen that at next to leading order in perturbation theory 
short time propagators begin to play a role. While this complication prevents the 
assignment of Wick contractions to orbit pairs of definite topology, the results obtained 
after integration over all temporal configurations remain universal (agree with the 
RMT prediction). We trust that the structures discussed above are exemplary for 
the behaviour of the ballistic <r-model at arbitrary orders of perturbation theory, 
i.e. that after integration over all intermediate times, each contraction contributing to 
(|3.16f) produces the universal result otherwise obtained by its zero dimensional analog 
equation (|3.15|) . 

4. Conclusions and Outlook 

In this paper, we have applied field theoretical methods to explore quantum 
interference corrections to the spectral form factor of individual chaotic systems. We 
have seen that the formation of the latter essentially relies on the fact that the ballistic 
cr-model — a field theory defined in classical phase space — is not capable of resolving 
structures on scales smaller than the Planck cell. This quantum uncertainty is an 
intrinsic feature of the model (viz. through the fact that the field degrees of freedom 
are multiplied by Moyal rather than by conventional products) and need not be added 

% While, in principle, the theory also permits the formation of short time propagators connecting 
two phase space points of a single vertex, these contributions are practically negligible: imagine a 
propagator Pt(x,x') returning after a short time to its point of departure (|x — x'| ~ ). Since t 
is much shorter than the Ehrenfest time, all other propagators departing from the concerned Hikami 
box will essentially follow the trajectory traced out by the return propagator, and, after a time t, 
also return to the departure region. In semiclassical language, we are dealing with an orbit that 
traverses a loop structure in phase space repeatedly. It is known, however, that for large time scales, 
the probability to find repetitive orbits is exponentially small (in the parameter exp(— Af)), i.e. short 
self-retracing contractions are negligible. 
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by hand as was done in previous approaches. In a manner that largely parallels 
the results of recent semiclassical analyses, the interplay of this uncertainty with the 
instabilities of the underlying classical chaotic dynamics leads to the formation of 
universal quantum interference corrections to the spectral form factor. 

The analysis above is perturbative in nature and, thus, limited to energy scales 
larger than the single particle level spacing. To advance into the perturbatively 
inaccessible regime u> < A (i.e. to prove the universality hypothesis in full) one 
would need to understand how the conspiracy of quantum uncertainty and classical 
instabilities damps out fluctuations inhomogeneous in phase space at time scales larger 
than the Ehrenfest time. The identification of a concrete mechanism effecting this 
reduction remains an open issue. 
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Appendix A. Regularization 

Throughout this appendix we use phase space coordinates x = (E,t,y), where E 
is the energy, t a time-like coordinate conjugate to energy and parameterizing the 
Hamiltonian flow through x, and y a (2/ — 2)-component vector of energy shell 
coordinates transverse to t. 

Appendix A.l. In-Shell Regularization 

As discussed in the text, the quadratic action of the model is controlled by the 
commutator [H, ] or, upon Wigner transformation, the series of operators [H, ] i— » 
ih{H, } + Y^n=i ^ 2 ™ +1 £ )(2 ™ +1) (5x), where D^ n+V > (<9 x ) is an operator of (2n + l)-th 
order in the phase space derivatives {d Xi }. When acting in a space of functions smooth 
on scales h, terms beyond the leading term (the Poisson bracket) are inessential and the 
quantum dynamics collapses to its semiclassical limit. Naively, one might hope that 
to achieve this reduction it suffices to choose the initial distributions in phase space 
sufficiently smooth. However, what complicates the problem is that the generator 
of classical evolution {H, } by itself leads to the dynamical buildup of singularities, 
no matter how smooth the initial distribution was. The point is that, due to the 
global hyperbolicity of the dynamics, we may locally identify truly expanding and 
contracting coordinate directions. Focusing attention on the latter, and linearizing 
the flow around a given reference trajectory, the equations of motion controlling the 
evolution of a phase space distribution p assume the form p = {H, p} = Xsd s p + . . ., 
where s is the coordinate that contracts strongest, A the corresponding Lyapunov 
exponent, and the ellipses indicate derivatives in other coordinate directions. After a 
time t ~ A -1 ln(5xo/h), where 5xq denotes the characteristic initial extension of the 
distribution, structures in the s-direction fluctuating on scales < h will have formed 
implying that derivatives D( 2n+1 ) acting in s-direction can no longer be neglected. 
One way to remove this complication [2H| is to add to the generator of classical time 
evolution an elliptic operator ~ Dd^, where D is constant. Indeed, it is straightforward 
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to show (by dimensional analysis or by explicit calculation) that for the regularized 
operator Xsd s + Dd 2 the initial contraction halts at a characteristic scale s ~ (D/X)i . 
Choosing D ~ ft 2 , it is guaranteed that the distribution will not build up structure 
on scales < ft, i.e. that the quantum corrections to classical dynamics are negligible. 
This motivates the addition of the regulatory contribution to the action. 

Appendix A. 2. Off-Shell Regularization 

In the main text (cf., e.g., equation l|3.12jl ). we have assumed that (a) all fields are 
defined on a shell of constant energy Eq and (b) the theory is local in the conjugate 
'time' coordinate t. To understand the meaning of this reduction, we need to recall 
the original definition of the field degrees of freedoms, Q(x), of the ballistic cr-model. 
According to Ref. 13 , Q = T^AT, where A = 6 E „( E o ~ H) <8> A and the 'delta 
function' 

projects on an energy window of width E av . (In the original reference, E av was 
identified with the energy window over which the two-point correlation function i|2.1[) 
is averaged, hence the subscript 'av'.) Integrating an action functional of these field 
degrees of freedom over all of phase space and absorbing Se„ into the integration 
measure, <ix i— > dx x Se^ v (Eq — H) ex (cbc), we obtain the 'energy shell' measure used 
in the text. To understand the energy dependence of the fields themselves, we write 
T(x) = 1 + W(x), where W(x) anti-commutes with A and, therefore, commutes 
with the function ^E av (^o — H). (Recall that all products are Moyal products, 
i.e. functions in phase space do not necessarily commute with each other.) Evaluating 
the latter condition, we obtain [W,S E ^(E - #)](x) = ih5' E ^(E - E)d t W + £>(ft 3 ). 
To rigorously fulfill this condition, we would need to require independence of W of 
the coordinate t along the flow through the phase space point x — obviously too 
strong a restriction. Instead, we will impose the weaker condition of approximate 
commutativity, | \[W, fe av (£ , — H)] | | 2 <C \\WdE av (E Q - H)\\ 2 , where the operator norm 
is defined as \\A\\ 2 = tr(A^A) = / dx|A(x)| 2 . It is straightforward to check — by 
explicit calculation or by dimensional reasoning — that the commutator is small in 
the above sense provided that d t W < (E av /h)W, i.e. the fields W have to be smooth 
in 'time direction' on scales ~ hjE aw .\ 

The above energy-time duality suggests the following interpretation of the theory: 
let us introduce a 'stroboscopic' picture of the particle dynamics wherein time scales 
smaller than a certain classical to need not be resolved. (E.g., in a billiard, to <C tf 
where tf is the time of flight through the system, etc.) All fields are smooth on 
scales to. This means that the width of the averaging window must be (at least) 
of order E av < h/to. In the classical limit, we indeed project onto a sharp 'energy 
shell'. (However, we do not know how to reconcile the condition of anti-commutativity 
with A with the condition of a mathematically sharp energy shell proposed in |17|.) 
Second, we require that the fields W(E) do not vary significantly over their narrow 
range [Eq — E av , E n +E av ] of definition (the so-called 'mode locking assumption' [T^]). 
This can be achieved by choosing the second order regulator derivative in E direction 

f Using equation 13.21 . it is also straightforward to show that field configurations W(E,t) which 
rigorously commute with 5^; av (^0 — H) do (a) vanish for energies outside a window of width _E av 
around and (b) have a bounded Fourier spectrum |e| < E^, where e is Fourier conjugate to the 
time variable r. 
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as ~ (ft/i ) 2 <9j;. We thus integrate over field configurations that are coarse grained 
over Planck cells of extension (to,E av = h/t ). As a result, the integral over the 
(E, t)-sector of the Moyal products can be carried out and we arrive at the theory 
independent of energies and local in time direction considered in the text. 



Appendix B. Perturbation Theory 

For completeness, we briefly summarize the contraction rules [13) employed in 
calculating integrals over products of B matrices: 

(tr(B(x)A) tr(BV)A')) B = ^-P w (x,x') tr(AA'), 



(tr(5(x)Ai?V)^)>B = ;^(x,x') tr(A) tr(A') 
(tr(B(x)A) tr(B(x / )A')) B = ^(x,x') tx{AA' T ), 
-{tr(B(x)AB(x')A')) B - ^(x,x') tv(AA' T ). 



(B.l) 



where A and A' are arbitrary fixed matrices. To compute the integral over an arbitrary 
product of traces of Z?-matrices, one first forms all possible total pairings B — B\ B — 
B, and B' — B\ and then computes individual pairings by means of fBljl . Each 
contraction reduces the number of matrices by two. Eventually, one obtains an 
expression ~ (tr 1)™ = (2i?)™- 2 (where all contributions with n > 2 vanish in the 
replica limit). 
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